Enhancement of entanglement transfer in a spin chain by phase shift control 
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We study the effect of a phase shift on the amount of transferrable two-spin entanglement in a spin 
chain. We consider a ferromagnetic Heisenberg/XY spin chain, both numerically and analytically, 
and two mechanisms to generate a phase shift, the Aharonov-Casher effect and the Dzyaloshinskii- 
Moriya interaction. In both cases, the maximum attainable entanglement is shown to be significantly 
enhanced, suggesting its potential usefulness in quantum information processing. 
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I. INTRODUCTION 

Transferring quantum information reliably and effi- 
ciently is an important task in quantum information pro- 
cessing. For example, quantum communication proto- 
cols, such as quantum key distribution (or quantum cryp- 
tography), usually require two (or more) distant parties 
to share entanglement of high quality to achieve tasks 
that are impossible in the regime of classical mechanics. 
Also in a typical situation we encounter in the standard 
quantum computation model, we need to couple two spa- 
tially separated qubits in order to perform two qubit uni- 
tary operations, e.g., a controlled-NOT gate. 

Most common approaches to this task include meth- 
ods with an information bus, guided ions (atoms), fly- 
ing photons, a sequence of swapping operations between 
neighboring qubits, etc. However, these methods require 
additional complexity in structures, manipulations and 
controls of the interaction between qubits, as well as re- 
peated conversions between the qubit state and another 
physical degree of freedom. Flying photons may be the 
best information carrier over macroscopic distances, but 
may not be so for microscopic scales of the order of, say, 
a few micrometers. This is why there has been intensive 
research activity in the past few years on quantum in- 
formation transfer via arrays/chains of stationary qubits 
that are interacting with their neighboring qubits. 

Typically, in previous studies of this topic, the quan- 
tum information channel consists of one or more chains 
of spin- 1/2 particles, each of which interacts with their 
nearest neighbors. The interaction between spins can be 
described by either the Heisenberg or the XY (or varia- 
tions of these, e.g., the XXZ) model with some relevant 
parameters for coupling strengths, anisotropy, etc. These 
types of models attract much attention because they are 
in principle sufficient for implementing quantum informa- 
tion processing: examples of proposed methods are those 
with quantum dots and particles trapped in an optical 
lattice. 



A number of important and interesting results have 
been reported in this research area: for example, sending 
quantum information through a spin chain without mod- 
ulation entanglement transport with an anisotropic 
XY model [2], perfect transfer by manipulating the cou- 
pling strengths [H, H|, near perfect transfer with uni- 
form couplings by a spatially varying magnetic field 
Fourier analysis-based quantum information encod- 
ing measurement-assisted transfer through two par- 
allel chains Q, specific realizations of the method in [l[ 
with superconducting qubit array [H, Q , perfect transfer 
by local measurements on individual spins [Iol |. and also 
transfer through a chain of coupled harmonic oscillators 

B 

The geometry of the spin chain could be more gen- 
eral in principle, but we will primarily consider a one- 
dimensional chain or ring of N spin-I/2 particles for sim- 
plicity. If we look at the chain as a whole, we can think 
of independent quantum states of collective modes, i.e., 
eigenstates of the Hamiltonian for the whole chain. The 
propagation of quantum information encoded in a spin 
can be thought of as the interference between all modes, 
which evolve independently. 

Hence, a naive strategy towards quantum information 
transfer of better quality would be to control the propa- 
gation of each mode, which will affect the (constructive) 
interference at a certain target site. Here we take this 
approach and consider the effect of changes in the en- 
ergy spectrum and dispersion relation of the mode (or 
'spin wave'), that induce a spin current in the chain. As 
a result of the induced change in the energy spectrum, 
there will also be a change in the time evolution of each 
mode and thus the interference between these modes. An 
advantage of this approach is that the pairwise coupling 
strengths Jij between i-th and j-th spins and the exter- 
nal electromagnetic field can be taken as a constant over 
the whole chain, unlike some schemes proposed before. 
That is, they do not have to be manipulated site by site, 
regardless of the starting and target sites for the transfer. 
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In this paper, we study the effect of a phase shift on 
the amount of transferrable entanglement through a spin 
chain. We shall primarily focus on the transfer of entan- 
glement, rather than the state itself, since entanglement 
is the key to achieving highly non-classical information 
processing. Besides, keeping the fidelity of the (entan- 
gled) two-spin system is harder than keeping the fidelity 
of a single spin. This can be stated more precisely as "the 
entanglement fidelity associated with a trace-preserving 
map is lower than or equal to the corresponding input- 
output fidelity of a subsystem" [12]. In short, attaining 
an entanglement transfer of high quality guarantees a 
state transfer of high quality. 

Further, we assume that the whole system is in the 
'one-magnon' state, in which the total number of up spins 
in the chain is one. This situation is simple enough to 
start our analysis with and is indeed reasonable when 
considering information transfer. Thus, all spins in the 
initial state, except the one that is the subject of the 
transfer operation, are initially down. As for the in- 
teraction between spins, we will consider the standard 
isotropic Heisenberg model. The Heisenberg model is ac- 
tually equivalent to the XY model under the one-magnon 
assumption, therefore our analysis will be applicable to a 
wide range of physical systems. Also, as assumed (mostly 
tacitly) in the literature listed above, we make another 
assumption not to let useful entanglement pass by the 
target site: the entanglement at the target site(s) can 
be extracted later on at will or can be retained for fur- 
ther operations, including entanglement distillation [13] . 
to maximize the efficiency of subsequent processes. 

We will consider two different ways of generating phase 
shifts. One is due to the Aharonov-Casher effect [Til ], 
while the other is induced by the Dzyaloshinskii-Moriya 
interaction [IH, [lj| • There are other means to generate 
a phase shift, or equivalently a (persistent ) sp in current, 
in a chain/ring, such as those reported in [171 [Til]. How- 
ever, we focus on the two above because these seem to 
be sufficient to demonstrate a significant entanglement- 
enhancing effect due to phase shifts. We discuss only 
chains of ring geometry because the effect is absent in 
open ended chains as shown in the Appendix. 



II. PAIRWISE ENTANGLEMENT IN SPIN 
CHAINS 

Let us start with a description of the entanglement 
between an arbitrary pair of spins in an ./V-spin chain 
within the one-magnon condition. The properties of the 
spin chain, such as its geometry, the nature of the inter- 
action, etc. could be any at this point. As there is only 
one up spin in total, the state of the whole chain at time 
t can be given by 



where 5+ = SJ + iSj is the raising operator defined with 
the spin-1/2 operators (a = x, y, z) for the j-th spin. 
The lowering operator is defined as its Hermite conju- 
gate, SJ — Sj — iSj. Throughout this paper, we will 
let |0) denote the spin-down state, while |1) is the spin- 
up state. As in standard entanglement transfer scenarios 
with the one-magnon assumption, the magnon is initially 
localized at a single site. With such initial conditions, we 
can identify the amplitudes ay(i) as the propagators, or 
the Green functions, from the point of view of wave me- 
chanics. They have all information to describe the time 
evolution of a magnon wave packet. 

In order to evaluate the pairwise entanglement in 
|^(t)), we employ the concurrence [l9j] as its measure. 
The concurrence C in a bipartite state p is defined as 

C := max{0, Ai - A 2 - A 3 - A 4 }, (2) 

where Ai are the square roots of the eigenvalues of ma- 
trix R in descending order. The matrix R is given as a 
product of p and its time-reversed state, namely 

R = p(a y ®a y )p*(a y ®a y ), (3) 

where a y is one of the standard Pauli matrices and the 
star denotes the complex conjugate. The concurrence 
takes its maximum value 1 for the maximal entanglement, 
while it is for all disentangled qubits. 

We will compute the concurrence Ci lt i 2 (t), between the 
li-th and /2-th spins, at time t by tracing out all spins 
except those two. Then, Ci lt i 2 {t) can be written as 

C lul2 (t) = 2\a h (t)\ ■ \a l2 (t)\. (4) 

Hence Eq. (0| states that the concurrence between the 
two sites in a chain can be expressed as a product of the 
absolute values of propagators. That the entanglement 
is determined by the propagators supports our approach 
towards enhancing the entanglement by phase shift, be- 
cause the propagators are naturally affected by change in 
dispersion relation caused by phase shifts. 

III. ENTANGLEMENT TRANSFER WITH THE 
AHARONOV-CASHER EFFECT 

First, we consider the Aharonov-Casher effect [lH as 
a physical mechanism that causes a phase shift in the 
collective modes. When a neutral particle with mag- 
netic moment p travels from r to (r + Ar) in the pres- 
ence of electric field E, the wave function of the particle 
acquires an extra phase, which is the Aharonov-Casher 
(AC) phase, 

A0 = — j / p, x E(x) ■ dx, (5) 
hc z jp 

in addition to the ordinary dynamical phase. The phys- 
ical origin of the AC effect is that a particle moving in 
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(a) 



(b) 




FIG. 1: Examples of configurations for the Aharonov-Casher 
effect in a ring-shaped spin chain. In (a), the z-axis is taken 
to be parallel to the direction perpendicular to the plane that 
contains the spin chain. With an electric field E, which is 
directed to the radial direction, the term /2 x E(x) ■ dx in 
Eq. ([5} takes its largest value. Alternatively, as in (b), the 
directions for the spin and the electric field can be swapped 
to have the same fl x E(x). The magnetic moment of spin 
eigenstate, | f) or | J.), is parallel to the radial direction. 



an electric field feels a magnetic field as well due to rela- 
tivistic effects: the AC effect is essentially equivalent to 
spin-orbit coupling. If there is no external field applied to 
the ring of spins, the dispersion relation should be sym- 
metric with respect to the zero wave number (fc = 0), 
i.e. Ek = E-k, due to the rotational symmetry of the 
system. However, if the accumulated AC phase along a 
ring does not vanish after a 2tt rotation, the dispersion 
relation will change as the applied field breaks the (spa- 
tial) symmetry. Consequently, the propagation speed of 
each mode will be affected and the concurrence between 
any two sites can be expected to change accordingly. 

Figure Q] sketches two possible configurations to have 
the AC phase effectively. The geometry in Fig. Ufa) i s 
very similar to that for electrons in an atom. An electric 
field diverges radially, and the z-axis is taken to be per- 
pendicular to the plane containing the ring. The term 
jl x E(x) ■ dx in Eq. ([5]) takes its largest value when 
a (quasi-) magnetic moment (magnon) travels along the 
chain. The electric field could be generated by, for ex- 
ample, a charge on a wire at the center of the ring. Al- 
ternatively, the directions for the spin and the electric 
field can be swapped as in Fig. |Tf b) to have the same 
/Z x E(x). 

We consider a Heisenberg chain of N spin-1/2 particles 
interacting ferromagnetically with their nearest neigh- 
bors. The phase acquisition due to the AC effect modifies 
the standard Heisenberg model Hamiltonian to [2(J, HH 
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shaped configuration is represented by periodic boundary 
conditions, i.e., N + 1 = 1. The Hamiltonian H can be 
diagonalized with the help of the Jordan- Wigner trans- 
formation [13, HH that maps spins under H to spinless 
fcrmions. The annihilation and creation operators for the 
fcrmion at site j are 



exp UiE^M j S i 



and c 



SJ exp 



j'-i 



(7) 



Under the one-magnon condition, the Hamiltonian H is 
now diagonalized as 



N r 

H = -E \ {?4<*» 

3 = 1 1 
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with a further linear transformation rjk = J^j 4'*kj c 3- Thfi 
energy eigenvalues E^ can be computed as 



E k =- cos(fc + 6) + 1 



N 



h 1 



N 
~2 



(9) 



where k = 2-xn/N with -N/2 < n < N/2, and <j) k j 
1 / ^/N e lk i . As the second and third terms of Eq. ([9]) are 
constant, we will omit them hereafter. A one-magnon 
eigenstate can be obtained accordingly with the form of 



vl as 



l*> :=^|0) 8Ar = 4rE elfe - ? ^ + l°) 



0N 



(10) 



The presence of the extra phase 8 is reflected only in the 
energy spectrum, while the expression for eigenstates is 
unchanged. The change in the dispersion relation by a 
phase shift is illustrated in Fig. [Ha). 

Note that applying the Jordan- Wigner transformation 
to the Hamiltonian © gives an additional term to Eq. 

®, -l/2{e w c\c N + e- ie cj v ci)(exp^^^ 1 c, t cO + 1), 
which is a result of the periodic boundary condition. We 
have already omitted this term in Eq. ([8]) since it equals 
zero as long as we consider the one-magnon state. This 
is because exp(i7r J^iLi c \ c i) + 1 = for any N. Also, 
the one-magnon condition makes the Heisenberg model 
equivalent to the XY model, as the interaction between 
(virtual) fermions after the Jordan- Wigner transforma- 
tion is absent under this condition in both models. 



where the interaction strength is taken as J = —1 for 
all neighboring pairs and h is the magnetic field, which 
is taken to be uniform and parallel to the z-direction. 
The phase change 9 between neighboring spins is given 



by Eq. © with 



and Ar 



rj+i 



The ring- 



A. Entanglement with an isolated spin 

Let us analyze the entanglement propagation along the 
spin chain. Suppose that at t = a physically isolated 
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(b) N=large 




k + 



FIG. 2: Energy spectrum Ek = — cos(fc+#) (in the units of J), 
(a) A non-zero 9 changes the dispersion relation represented 
with open circles (no phase shift) to the one with filled circles 
(0/0). The number of sites N is taken to be 8. (b) The 
energy spectrum when TV is large. 



(the 0th) spin and the first spin are maximally entangled 
as (|01) + \10))/y/2 and the rest of the spins in the chain 
are all in |0). Thus, the initial state of the whole system 
can be expressed as 



l*(0)> = 



v/2 



|0)o 



lk \k) 



|l)o|0) 



(11) 



Hence, we find the state at time t, taking h = 1, as 



l*M> - T2 



lyjexplitO-lJ-iEjt) 



|l)o|0) 



®N1 



(12) 



Figure [3] depicts the process we consider: Fig. [3ja) shows 
the initial correlation in |^(0)), and Fig. [3^b) is the 
desired goal of our entanglement transfer operation. 

Now we can evaluate the entanglement between the 
0th and Z-th spins. Equation (|12[) can be written in the 
form of Eq. |T]) with amplitudes 



/(*) = 



1 



V2N^ 



exp [ik(j - 1) - iE k t] . 



(13) 



Because the 0th spin is not interacting with other spins, 
we can take |an(i)| = l/\/2 for all t. From Eq. (0| we 



(a) initial state 

entanglement 



(b) desired final state 
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FIG. 3: Schematic picture of entanglement transfer in a spin 
chain, (a) The configuration of the Heisenberg ring and the 
entanglement it has at t = 0, with a spatially separated sys- 
tem that is here represented as the 0th spin, (b) The ideal 
goal of the entanglement transfer operation: we wish to trans- 
fer as much entanglement with the 0th spin as possible to a 
specific (target) spin in the chain. 



obtain Cq\ (t) — y/2\oti(t)\, where the superscripts on C 
denote the initially entangled pair, fn the limit of large 
N, this takes a simple analytical form 



e -<('-i)(«-f >./,_!(*) =|Ji_i(f)|, (14) 



with the Bessel function of the first kind J u (x). The 
effect of the AC phase 9 disappears in this limit, since 
the energy spectrum becomes continuous and displacing 
all modes by 9 does not change the overall dispersion 
relation. This can be clearly seen in the plot in Fig. 
[2jb). In other words, 9 appears only as a common phase 
factor for all modes, e - ^' -1 ^ 9- ^, thus there is no 9- 
dependence in \atj\. 

An example of the plots of concurrence as a function 
of t and 9 is shown in Fig. QJ which is the plot of C ' 3 (i) 
for N = 5. The phase 9 is the same as that in Eq. 
(IH), that is, the phase a magnon acquires when hopping 
from the j-th to the (j + l)-th site. Some improvement 
in the amount of transferred entanglement due to the 
nonzero phase shift is evident: the maximum concurrence 
when 9 = is C^=° = 0.647 (at t = 59.05), and when 
9^0, C max can reach as high as 0.996 (at t — 23.71 
and tan# = 1.376) in the region we have calculated, i.e., 
t E [0,200] and 6 £ [-tt.tt]. 

Figure [5] shows the comparison between the maximum 
values of concurrence C max with and without phase shift 
for various (N, I) from (3,2) to (13, 13), where N is the 
total number of sites and I is the site where the concur- 
rence is evaluated. Plotted are the highest values found 
numerically in the range of < t < 200. The horizon- 
tal axis represents (N,l). The blue lines with diamond 
markers are for the geometry with an isolated spin, and 
the red lines are for the entanglement transfer when the 
initial entanglement is held by a pair in the chain. The 
latter case will be discussed in the following subsection. 
For both cases, the open markers show the maximum 
concurrence CT 1 )' 1712 when 9 — 0, while the filled ones 
mark maxjC™ 1 ^™ 2 } in the range of — it < 9 < tt. 
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FIG. 4: An example of plots of the concurrence C transferred 
in a five-spin chain (ring) in the presence of a phase shift. The 
initial entanglement (of the form (|01) + |10))/V2) is in the 
first spin of the ring and an isolated (0th) spin. The phase 
shift 9 is an extra phase a magnon acquires when travelling 
to a neighboring site. The concurrence C is evaluated for the 
pair of the 0th and third spins. Units for t and 6 are h seconds 
and radian, respectively, where h is the Planck constant. 



The enhancement of entanglement by the phase shift 
is clearly seen in Fig. [5l The transferred entanglement is 
significantly increased by nonzero phase shifts for many 
values of (N,l). Yet, the degree of enhancement varies 
because quite an effective constructive interference can 
occur even when 9 = for some (N, I). 



B. Entanglement in a pair of spins in the chain 

Let us now look at an alternative scenario of entangle- 
ment transfer. Instead of being entangled with an iso- 
lated (0th) spin, both spins of an entangled pair can be 
in the same chain. Considering the effect of the phase 
shift mentioned above, we can naturally expect some im- 
provement in the efficiency of transfer in this scenario as 
well. We already have the expression Eq. Q for the con- 
currence at time t between the Z-th and (I + l)-th spins 



ol mi, in 2 



(01 



! (t)|, where mi and 



m.2 denote the initial sites that are entangled. Factors 
a™ 1 '™ 2 ^) are given explicitly by 

«r m2 W = ^E(exp[ifc(3-m 1 )-i£?*t] 

+ exp[ik(j - m 2 ) - iEkt]) . (15) 

The comparison of the maximum concurrence between 
the cases with and without the phase shift is shown in 
Fig. [5] with red triangular markers. For the plot in Fig. 
[5J the initial entanglement of the form of (|01) + |10))/\/2 
is assumed to be in the first and second spins and the rest 
are in |0). 



As N tends to infinity, C™V™ 2 (t) approaches the form 

C™r 2 W = \jh- mi (t) + e- l ^- m ^ e -"^J h ^ m2 (t) 
x J h - mi (t) + e-^- m ^ 6 -*IVj l2 ^ m2 {t)\ . (16) 

Interestingly, unlike Eq. (|14p. there is still a dependence 
on 8, regardless of Vs and m's. Nonzero 8 can indeed 
always increase the maximum attainable pairwise entan- 
glement in a long chain, no matter which pair is initially 
entangled, and no matter which pair we evaluate the con- 
currence for. 

We can see in Fig. [S] that the degree of enhancement 
is larger when the two initially entangled spins are in 
the chain, compared with the case of entanglement with 
an isolated one. This difference can be understood intu- 
itively when N is large: the physical reasoning for small 
N's is essentially the same. If we look at a propagator 
from a single site, the phase shift 8 gives an almost com- 
mon displacement to the phase of all modes at any single 
site as we have seen in Eq. (fT^)l . When the entangled 
pair, mi and 7772, is embedded in the chain initially, there 
are two independent propagators stemming from the two 
sites. Since these propagators have different phase dis- 
placements at the same site, say the Z-th, the interference 
due to non-zero 8 still occurs. As a result, there remains 
a dependence of the concurrence on the phase shift. 



IV. SPIN SHIFT INDUCED BY THE 
DZYALOSHINSKII-MORIYA INTERACTION 

Phase shifts can be generated by a different type of 
interaction, that is, the antisymmetric exchange interac- 
tion, J a.k.a. the Dzyaloshinskii-Moriya (DM) interaction 
[HI, EH, in solids. The DM interaction could be quite 
significant in some solid-state-based qubit systems, such 
as quantum dots (24|. 

The Hamiltonian for the DM interaction between two 
spins, say 1 and 2, can be written as 



H 



DM 



d-(S 1 



(17) 



where d is the coupling vector that reflects the anisotropy 
of the system. Assuming that only the z-component of 
d has a nonzero value, i.e., d z ^ 0,d x = d y = 0, and 
that all components are constant along the chain, we can 
write the total interaction Hamiltonian and its spectrum 
as (omitting the terms that give only a constant bias) 



H 



3 
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— 'J2cos(k + <f>)r)lri k , (18) 
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I 2 2 2 2 2 5 2 5 25 258258 258 258 11 

FIG. 5: Comparison of the maximum concurrence attainable C ma x between the cases with/without phase shift. The blue 
plots (diamonds) show the concurrence between the 0th and Z-th spins when the 0th and first spins have the entanglement 
(|01) + |10))/V2 at t = 0. The filled and open diamonds correspond to nonzero and zero phase shift, respectively. The red plots 
(triangles) are for the concurrence between the Z-th and (Z + l)-th spins when the first and the second spins are initially entangled 
in the same form. The horizontal axis represents the total number N of sites in the chain and the location I (2 < I < N) at 
which the concurrence is evaluated. 



where </> = tan -1 d z . Thus, the difference from the energy 
spectrum under the AC effect, Eq. ©, is only the factor 
1 / cos </> for the energy eigenvalues. As this factor is inde- 
pendent of fc, it only rescales the energy spectrum linearly 
(for a given cj)) , hence increases the speed of propagation 
of all modes by 1/ cos0. Consequently, the maximum 
concurrence attainable stays the same as that in the AC 
effect case, though the time at which the maximum is 
achieved should be rescaled as well. All quantitative re- 
sults in the previous section are valid for this system with 
H dm if t is replaced with tj coscj). 



V. SUMMARY AND OUTLOOK 

We have investigated the effect of externally induced 
phase shifts on the amount of entanglement that can be 
transported in a spin chain. As we have qualitatively an- 
ticipated in the Introduction, these phase shifts can sig- 
nificantly enhance the efficiency of entanglement trans- 
fer. Although we have only studied two shift generating 
mechanisms, we believe that phase shifts are useful in 
quantum information processing, particularly for short- 
distance transfers, regardless of the mechanism. Also, 
we have found that there is an interesting clear differ- 
ence in the response to nonzero shift when the chain is 
sufficiently long. 

In the AC-effect-related experiment, there could be 
a difficulty in providing an electric field that is intense 



enough. A rough calculation gives an estimate of the nec- 
essary strength of the electric field of at least 10 7 V/m to 
have a meaningful phase shift, if the system size is of the 
order of a /zm. Nevertheless, such a strong electric field 
can be realized by two dimensional electron gases formed 
in heterostructured SiGe, GaAs, or other types of III-V 
materials. Furthermore, the so-called band-enhancement 
of spin-orbit coupling in crystals [25| could be useful to 
have a substantial AC effect. 

Despite a number of technical difficulties, some qubit 
arrays, in which the effect of the phase shift can be 
observed, could be fabricated with present-day technol- 
ogy. For example, consider an array formed with charge 
qubits, one type of superconducting qubits [2(|. Quan- 
tum information in a charge qubit is represented by the 
number of excess Cooper pairs in the superconducting 
Cooper-pair box. When neighboring qubits are coupled 
via a Josephson junction, the effective interaction can 
be described by the XY model [27j. As the informa- 
tion carrier in this case is a pair of electrons, a phase 
shift can be induced to the wave function of the pair by 
the Aharonov-Bohm effect. A magnetic flux <E>b thread- 
ing through a ring formed by charge qubits with the 
Josephson-junction-mediated coupling would generate a 
phase shift (e/hc)§B for the wave function of the Cooper 
pair. Then we could expect the same effect discussed in 
this paper. Nonetheless, observing this effect is by no 
means straightforward: all technical problems from nano- 
structure fabrication to measurement method should be 
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addressed. We shall leave these challenging experimental 
problems for future investigation. 
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Appendix: Note on open ended chains 

As a geometry for information transfer, open ended 
linear chains may look more natural. If the dispersion 
relation can be affected by the phase shift in the case 
of open ended chain, then the amount of entanglement 
transferred can be expected to change as well. However, 
this is not the case. Let us briefly look at this. 

Suppose that a chain of N spins is placed in a uniform 
electric field as in Section IIIII The Hamiltonian is the 



same as Eq. ([6]), but instead of the periodic boundary 
condition we have the open boundary condition (OBC), 
c*o = af/+i = 0. The one-magnon eigenstates are given 
by 

where k = nn/(N + 1) with -N/2 < n < N/2. Corre- 
sponding energy eigenvalues are 

^ BC = -cosfc. (20) 

Clearly, the phase shift 9 has no effect on the energy 
spectrum and thus the propagation speed of each mode. 
Therefore 9 causes no change in the concurrence between 
any two sites. 

This result can also be paraphrased in the following 
way. The phase shifts at all sites can be cancelled by a 
product of local gauge transformations, IL- exp[ij9(Sj + 
1/2)], in the case of open ended chains. The same can- 
cellation cannot be made for ring-shaped chains because 
of the accumulated phase along the chain. 
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